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ABSTRACT 

In this paper, we present the Fermi All-sky Variability Analysis (FAVA), a tool to systematically study the variability 
of the gamma-ray sky measured by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. 

For each direction on the sky, FAVA compares the number of gamma-rays observed in a given time window to 
the number of gamma-rays expected for the average emission detected from that direction. This method is used 
in weekly time intervals to derive a list of 215 flaring gamma-ray sources. We proceed to discuss the 27 sources 
found at Galactic latitudes smaller than 10° and show that, despite their low latitudes, most of them are likely of 
extragalactic origin. 

Key words: binaries: general - BL Lacertae objects: general - catalogs - galaxies: active - stars: flare - surveys 
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1. INTRODUCTION 

In 1844, the astronomer F. W. A. Argelander performed one 
of the first systematic studies of the variability of the night 
sky. He laid the study of variable sources “most pressingly 
on the heart of all lovers of the starry heavens, to perform an 
important part toward the increase of human knowledge, and 
help to investigate the eternal laws which announce in endless 
distance the Almighty power and wisdom of the Creator” (Percy 
2007). Nowadays astronomers are not as poetic, but time has 
provided us with exceptional instruments for the quest. 

In this paper, we present a systematic study of the temporal 
variations of the gamma-ray sky measured by the Large Area 
Telescope (LAT) on board the Fermi satellite. The gamma- 
ray sky above 100 MeV is dominated by the Galactic diffuse 
emission, which originates from cosmic-ray interactions with 
interstellar matter and photon fields (e.g., Ackermann et al. 
2012b). Additionally, an isotropic diffuse gamma-ray emission 
is detected, and it is the strongest source of emission at high 
Galactic latitudes (Abdo et al. 2010d). Both diffuse components 
are expected to be stable over the duration of the Fermi mission. 
On top of this background, 1873 gamma-ray sources have been 
detected during the first two years of the Fermi mission and 
reported in the second Fermi - LAT catalog (2FGL; Nolan et al. 
2012). Out of these sources ~24% are found to be variable on 
monthly timescales. The vast majority of the variable sources 
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are associated with active galactic nuclei (AGNs), which are 
known to be variable across the electromagnetic spectrum (e.g., 
Aharonian et al. 2009; Ackermann et al. 201 la, 201 lb). 

Variability in gamma-rays has so far been established only 
for a few sources in our Galaxy. Orbital modulation and isolated 
flares have been reported from seven X-ray binaries, in which 
a neutron star or a black hole orbits a massive companion 
(Aharonian et al. 2005, 2006; Albert et al. 2006; Abdo et al. 
2009b; Hinton et al. 2009; Sabatini et al. 2010; Ackermann 
et al. 2012c). Variable gamma-ray emission has also been 
reported from the direction of 17 Car, a massive star Wolf-Rayet 
binary (Tavani et al. 2009; Abdo et al. 2010b; Reitberger et al. 
2012). Recently, flaring gamma-ray emission has been found 
for two new source classes: Nova explosions (Abdo et al. 
2010c) and the Crab pulsar wind nebula (Abdo et al. 2011c; 
Tavani et al. 2011; Buehler et al. 2012). The latter was thought 
to be a stable gamma-ray emitter, but was discovered to be 
flaring with the method described in this paper. In addition, 
several other gamma-ray transients have been detected near 
the Galactic equator, but they are likely associated with distant 
AGNs (Vandenbroucke et al. 2010; Cheung et al. 2012a). 

The reasons for the detection of so few variable gamma- 
ray sources within our Galaxy remain unclear, whether due 
to astrophysical reasons, due to the statistically limited flux 
sensitivity, or due to systematic difficulties of detecting them 
owing to uncertainties in the modeling of the strong foreground 
of the Galactic diffuse emission. In this paper, we present a new 
method developed to search for transients in the gamma-ray sky 
that does not require a diffuse emission model. We first describe 
the methods and proceed to assemble a list of flaring gamma-ray 
sources seen over the sky during the first 47 months of the Fermi 
mission. We then focus on the sources detected at low Galactic 
latitudes, as they may be of Galactic origin. 
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2. THE FERMI ALL-SKY VARIABILITY ANALYSIS 

Flux variability of LAT sources is usually studied with a 
maximum likelihood analysis, in which parameters of a model 
describing the point sources and diffuse gamma-ray emission in 
a given region of the sky are jointly optimized. The sensitivity of 
this approach is often limited by the uncertainties of the diffuse 
emission modeling, particularly in the Galactic plane (Nolan 
et al. 2012). Small inaccuracies in the instrument response 
functions can lead to time-dependent residuals which depend 
on varying observation conditions, e.g., off-axis angle of the 
sources or orbital position of the LAT, limiting variability studies 
(Ackermann et al. 2012a). A further limiting factor of the 
likelihood approach is that it is computationally intensive; it 
is currently difficult to perform variability studies in different 
time and energy windows over the entire sky. We therefore 
developed the Fermi All-sky Variability Analysis (FAVA), in 
which we search over a grid of regions on the sky for deviations 
from the expected flux based on the long-term average. While 
this approach is less sensitive than a likelihood analysis, it has 
three main advantages. 

1 . The analysis is independent of any model for the diffuse 
gamma-ray emission. The diffuse emission is expected to 
be constant over the time of the Fermi mission. It therefore 
cancels out in the comparisons between the number of 
expected and measured gamma-ray events. 

2. The analysis is computationally inexpensive, allowing us 
to blindly search for flux variations over the entire sky. 
The analysis is therefore unbiased, treating every direc- 
tion on the sky equally, potentially yielding unexpected 
discoveries. 

3. No assumptions are made about the spectral shapes of the 
gamma-ray sources. Negative flux variations are treated the 
same way as positive ones. (Throughout this paper, we refer 
to both positive and negative variations from the mean as 
flares.) 

We applied FAVA to the first 47 months of Fermi observations 
(2008 August 4 to 2012 July 16 UTC), in weekly time intervals. 
The total number of weeks is 206. We considered two ranges 
of gamma-ray energy, E > 100 MeV and E > 800 MeV, 
to increase the sensitivity for spectrally soft and hard flares, 
respectively. We used the P7SOURCE_V6 event selection and 
only considered events with a zenith angle smaller than 95°, to 
limit contamination from the gamma-ray emission of the Earth 
atmosphere, which is time variable in sky coordinates. 

We generate measured and expected count maps with a resolu- 
tion of 0.25 deg 2 per pixel. The maps are smoothed by assigning 
to each pixel all events that were detected within a distance cor- 
responding to the 68% containment radius of the point-spread 
function (PSF). The pixel positions are characterized in spheri- 
cal coordinates by 0 and 9. The number IV exp (0, 9) of expected 
events in 1 pixel is derived from the number !V tot (0, 9) of events 
observed from the same direction over the first 47 months of 
observations. As the PSF depends on the photon energy E and 
on the incidence angle a with respect to the LAT (Ackermann 
et al. 2012a), we integrate over these parameters: 


A exp (0, 9) 


E E x 


<f k 0M) 
€%{</>, 0 ) ’ 


( 1 ) 


are a statistical fluctuation of the expected value, based on 
Poisson statistics. These probabilities are then converted to 
significance in units of Gaussian sigmas for easier visualization. 
The procedure is illustrated in Figure 1 for one example week 
in 2009 February, during which the first flare from the Crab 
Nebula was seen in LAT data. 

The sensitivity of FAVA to detect flares varies with the po- 
sition in the sky, due to the anisotropy of the diffuse emis- 
sion backgrounds. Additionally, the sensitivity depends on the 
energy spectrum of the flaring gamma-ray source. Typically, 
gamma-ray sources have photon indices between 1.5 and 2.5 
in the Fermi-LAJ energy range (Nolan et al. 2012). The sen- 
sitivity for both of these cases is shown in Figure 2. Flares 
with a photon index greater than ~2 are typically detected 
at higher significance in the low-energy maps, while those 
with a smaller index are detected more significantly in the 
high-energy maps. 

The accuracy of FAVA was tested on simulations of a con- 
stant sky for 36 months of observations. The simulated sky 
was composed of the Galactic (gal_2yearp7v6_v0.fits) and 
isotropic (iso_p7v6source.txt) diffuse emission models 64 and 
point sources. The latter were generated with random coordi- 
nates on the sky and according to the flux distribution derived 
in Abdo et al. (2010e). The distribution of the significance a of 
flux variations is displayed in Figure 3 for the simulations and 
the real data. In the low-energy interval, on average 154 events 
are recorded per sky pixel in each week. The simulated signifi- 
cance distribution is therefore expected to be close to Gaussian 
in the low-energy band. Indeed, a fit to the simulations shows 
that a Gaussian model with a standard deviation of 1.00 and 
a mean of 0.06 describes the distribution well, as shown in 
Figure 3. In the high-energy band, an average of 3.6 events 
are recorded per sky pixel. The simulated distribution of flux 
variations is therefore not expected to be Gaussian. In partic- 
ular, typically only a few counts are detected in several pixels 
at higher Galactic latitudes, leading to a small-scale structure 
around \o\ < 1 in the significance distribution. A disagreement 
is visible in this region between the data and simulations. How- 
ever, for |oj >2 the qualitative agreement with expectations 
is good also in the high-energy band, we therefore expect no 
significant biases in the detection of flares at high significance. 
As at low energies, one can see that the difference between 
simulated and real significance distributions increases with in- 
creasing cr, corresponding to the real flux variations emerging 
above the statistical background. 

The good agreement of the statistical fluctuations in the low- 
energy band of the simulated constant sky with expectations 
shows that systematic effects are small in the simulations. 
However, as seen for the discrepancy in the |cr| < 1 interval 
in the high-energy band, additional effects might be present in 
the real data, as the simulations do not take into account possible 
sources of systematic errors. In particular, they neither account 
for any disagreement between the simulated and real instrument 
responses, nor for varying background levels due to residual 
cosmic rays mistakenly classified as gamma-rays. It is difficult 
to assess these systematic effects from the data, as, e.g., small 
flux variations might be present throughout the gamma-ray sky 
due to variable background sources. However, we can set upper 
limits on possible systematic errors by looking at presumably 
constant sources such as pulsars. Analysis of the brightest 
ones, the Vela and Geminga pulsars, shows that their relative 


where e week and e tot are weekly and total exposures. We 
proceed to calculate the probability that the observed counts 
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Figure 1. Illustration of the map generation in FAVA. The week shown covers the time interval MJD 54864.655-54871.655 (2009 February 2-9). The measured 
counts map is shown in the first row, for energies >100 MeV (left) and energies >800 MeV (right, white color indicates that no counts were detected from this area). 
The middle panels show the expected counts from the average emission observed during the first 47 months of Fermi observations. The generation of measured and 
expected counts is explained in the text. The third row shows the significance of the flux variations. As an example, the position of the Crab Nebula is indicated by a 
star in the lower left panel. Its flux was increased compared to average with a significance >5 a for energies > 100 MeV during this week. The flare is only detected in 
the low energy range, as the energy spectrum was very soft (photon index ~ 3.5; see Abdo et al. 201 lc). An example of a flare detection with a negative flux variation 
is given by the blazar 3C 454.3. Its flux was lower than average for both the > 100 MeV and >800 MeV energy ranges with a significance >8a. Figures are shown in 
Galactic coordinates in a Hammer-Aitoff projection. Note that the color scales were adjusted to different ranges for the low- and high-energy bands. 



Figure 2. Minimum flux increase required for a flare detection in FAVA in a one- week interval with a significance of 5.5a in the low-energy or the high-energy 
band. Detection above these thresholds results in inclusion in the list of flaring sources, as discussed in Section 3. The figures are shown in Galactic coordinates in 
a Hammer-Aitoff projection. The flaring source is assumed to have a power-law spectrum in energy. The left panel shows the sensitivity for spectrally hard flares 
(photon index 1.5) and the right panel shows the sensitivity for soft flares (photon index 2.5). 


count variations are compatible with a steady flux within <5% 3. LIST OF FLARING GAMMA-RAY SOURCES 

on weekly timescales. Systematic errors of FAVA for relative 
flux variations of bright sources are therefore comparable to 
those for the standard Fermi-LXT analysis (Ackermann et al. 

2012a). 
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After calculating the significance maps for all weeks, we scan 
them for significant flares. To determine their positions, we use 
the peak finding algorithm described in Morhac et al. (2000). 
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Significance Significance 

Figure 3. Distribution of the significance of flux variations for the low-energy (left) and high-energy (right) intervals. The integral of both distributions has been 
normalized to one. Open circles show the simulations of a constant gamma-ray sky. Crosses show the measurement of the first 47 months of Fermi observations. The 
solid line in the left panel shows the best-fit Gaussian model for the simulated distribution (see the text). 


Based on these detections, we build a list of flaring sources. 
For this we only consider flares with significances greater than 
5.5er in the low-energy or high-energy band. This threshold 
was chosen so that the number of false flare detections due to 
statistical fluctuations is expected to be ~1 over the 206 weeks 
that were analyzed/ 1 ' 1 No flares were detected in the simulations 
of three years of data for a steady sky above this threshold. The 
number of false flares expected for 47 months of data is therefore 
<3 at 90% confidence. Additionally, we only consider flares that 
occurred far away from the average position of the Sun in the 
corresponding week. The Sun is a bright gamma-ray source that 
moves along the ecliptic by ~ 1 ° per week (Abdo et al. 201 lb). 
We therefore only considered flares at a distance from the Sun 
>12° and >8° in the low- and high-energy bands, respectively. 
Finally, we merge low- and high-energy flares detected in the 
same week, if they are coincident in position within 3°, relating 
them to the same flaring event. For the position of the latter we 
use the position of the high-energy detection, due to its higher 
accuracy, as will be discussed in the next paragraph. A total of 
1419 flares that fulfill the mentioned criteria were detected. Out 
of these, 645 and 175 are detected in the low- and high-energy 
interval only, respectively. The remaining 599 flares are detected 
simultaneously in both energy bands. 

To estimate the position accuracy achieved by the peak finding 
algorithm, we analyzed the distribution of flares around known 
flaring gamma-ray sources. As a source sample we chose the 
249 sources flagged as most variable in 2FGL (a variability 
index >83.2 in Nolan et al. 2012). The resulting distribution 
of flares per solid angle is shown in Figure 4. We assumed 
that the reconstructed position of the peak finder follows a 
Gaussian distribution plus a constant background term from 
flares not associated with the sources in this representation. The 
best-fit model shown in Figure 4 represents the data in good 
approximation. We proceeded to calculate the distances within 


65 The number of trials can be approximately estimated as the total sky area 
divided by the area of the PSF. Above 100 MeV the 68% containment radius of 
the event-averaged PSF of the l enm-\ .AT is ~3°. We therefore have 
~41253/(7r x 3 2 ) x 206 = 300558 independent tests in the sky. This results in 
~0.01 expected false positives above 5.5<r. For the high-energy maps the 
average PSF is ~0?6, resulting in ~0.3 expected false positives above the same 
threshold. 



Figure 4. Distribution of angular distances between flares exceeding the 
significance criteria for source detection and known variable gamma-ray 
sources. The solid and dashed lines show the best-fit parameterizations for 
flares with a detection in the high-energy band, and only in the low-energy 
band, respectively (see the text). 


which 68% of the flares are contained in this parameterization. 
At low energies the radius is 07 8 and at high energies it is 0:6. 
We verified that no systematic offset is present in the position 
estimation in any coordinate direction. 

Most variable gamma-ray sources, such as AGNs or X-ray 
binaries, are known to have recurring flares. We therefore group 
the detected flares, and associate closely located flares to a 
single common flaring source. For this purpose we used a 
Minimum Spanning Tree (MST; Nea et al. 2001). We first group 
the flares detected at high energies because their positions are 
better determined. We build the MST for these flares and merge 
neighboring flares with a distance of less than 2° in the spanning 
tree. We proceed to associate flares detected only at low energies 
to the ones found at high energy if their distance is less than 3°. 
Finally, we build the MST of the low-energy flares that were not 
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Figure 5. FAVA sources are shown in Galactic coordinates and a Hammer-Aitoff projection. Red crosses mark sources for which at least one flare was detected in the 
high-energy band. Sources that were detected only at low energies are marked by a yellow X. The colored image in the background shows the maximum significance 
<7 m ax detected in each pixel either in the >100 MeV or >800 MeV energy bands during the first 47 months of Fermi observations. The lower three panels show the 
region of Galactic latitude within 10° of the equator, the region enclosed by the dashed lines in the upper panel. 


associated and merge neighboring flares with a distance of less 
than 3° in the MST. 

The position of a FAVA source is found by averaging over the 
positions of its flares. If the source is detected at high energies 
we use only the positions of these flares due to their better 
positional accuracy, otherwise the low-energy flares are used. 
The position error of the source is obtained by propagating 
the positional errors of the flares included. In addition to this 
statistical error, there is a systematic error that can arise from 
false associations of flares to a source, as well as the finite 
binning of the sky maps. We estimated this error to be smaller 
than r sys = 0. 1 by comparing the position of the FAVA sources 
to those in the list of variable 2FGL sources used previously. 
We assume r sys to be the systematic error on the source 
positions. 

A total of 215 sources are detected by FAVA. Out of these, 33 
are detected at low energies only. Flares related to negative flux 


variations from the average emission are found for 22 sources, 
often during periods of quiescent emission. All of the latter 
also showed positive flares. No source was found which flared 
only due to a negative flux variation. Each FAVA source is 
referred to by its identification number composed of the right 
ascension in hours and minutes and the declination in degrees 
of the source (1FAV HHMM-DD). The positions of all sources 
in the sky are displayed in Figure 5. We produced light curves 
of relative flux variations with FAVA for all sources and made 
them publicly available online. 6 ' 1 One example light curve is 
shown for the position of the high-mass X-ray binary Cyg X-3 in 
Figure 6. 

We looked for associations of FAVA sources with previously 
known variable LAT sources. We searched for counterparts 
within radius R s , which is defined as the 99% statistical error on 


66 https://www-glast.stanford.edu/pub_data/585 
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Table 1 

List of FAVA Sources 


Ackermann et al. 


1FAV ID 

1 

(°) 

b 

(°) 

^68 

(°) 

N 1 

n L 

N f 

neg 

Rs 

(°) 

LAT Assoc. 

ATel 

Assoc. 

J0019-05 

102.0 

-66.9 

0.8 

1 

0 

0 

1.8 

2FGL J0017.6— 0510 


PMN J0017— 0512 

J0029-55 

309.7 

-61.8 

0.8 

1 

0 

0 

1.8 

2FGL J0032.7— 5521 



J0031-02 

111.1 

-64.7 

0.6 

1 

1 

0 

1.4 




JO 104+5 8 

124.6 

-4.5 

0.3 

7 

5 

0 

0.7 

2FGL J0102. 7+5827 


TXS 0059+581 

JO 112+23 

129.1 

-39.6 

0.6 

1 

1 

0 

1.4 

2FGLJ01 12. 1+2245 


S2 0109+22 

J01 12+61 

125.5 

-1.5 

0.3 

7 

3 

0 

0.8 

2FGL J0109. 9+6132 


TXS 0106+612 

J01 13+32 

128.4 

-30.4 

0.4 

4 

2 

0 

1.0 

2FGLJ01 12.8+3208 


4C 31.03 

J0115-11 

144.2 

-73.4 

0.3 

5 

5 

0 

0.7 

2FGL JO 11 6.0— 1134 


PKS 01 13-118 

JO 139+47 

131.2 

-14.6 

0.8 

1 

0 

0 

1.8 

2FGL J0136. 9+4751 


OC 457 

J0203-17 

186.2 

-70.6 

0.8 

1 

0 

0 

1.8 

2FGL J0205.3— 1657 


PKS 0202-17 

J0203+15 

147.5 

-43.9 

0.6 

1 

1 

0 

1.4 

2FGL J0205. 0+1514 


4C +15.05 

J0203+30 

140.9 

-30.0 

0.3 

5 

5 

0 

0.7 




J0209-51 

276.6 

-61.8 

0.3 

10 

5 

0 

0.7 

2FGL J0210.7— 5102 


PKS 0208-512 

J0212+10 

153.6 

-47.9 

0.3 

5 

5 

0 

0.7 




J02 15+01 

161.2 

-54.7 

0.6 

1 

1 

0 

1.4 

2FGL J0217. 9+0143 


PKS 0215+015 

J0219-27 

220.6 

-70.4 

0.6 

2 

1 

0 

1.4 




J02 19+35 

142.3 

-23.8 

0.6 

1 

1 

0 

1.4 

2FGLJ022 1.0+3555 


S4 0218+35 

J0225+42 

140.6 

-16.7 

0.3 

4 

3 

0 

0.8 

2FGL J0222. 6+4302 


3C 66A 

J0227-56 

279.1 

-56.0 

0.6 

1 

1 

0 

1.4 




J0228-36 

243.5 

-67.4 

0.6 

3 

1 

0 

1.4 

2FGL J0229.3— 3644 


PKS 0227-369 

J0234-61 

283.5 

-51.4 

0.3 

3 

3 

0 

0.8 

2FGL J0237. 1—6136 


PKS 0235-618 

J0237+60 

135.5 

0.4 

0.4 

4 

0 

2 

0.9 

2FGL J0240.5+61 13 


LS 1+61 303 

J0238+16 

156.9 

-39.2 

0.1 

16 

14 

1 

0.4 

2FGL J0238. 7+1637 


AO 0235+164 

J0238+28 

149.6 

-28.7 

0.2 

13 

7 

0 

0.6 

2FGL J0237. 8+2846 


4C +28.07 

J0246-46 

261.4 

-60.1 

0.4 

3 

2 

0 

1.0 

2FGL J0245.9— 4652 


PKS 0244-470 

J0251-22 

210.2 

-62.5 

0.3 

7 

4 

0 

0.7 

2FGL J0252.7— 2218 


PKS 0250-225 

J0302-24 

215.2 

-60.4 

0.3 

3 

3 

0 

0.8 

2FGL J0303.4— 2407 


PKS 0301-243 

J0308+10 

169.2 

-39.8 

0.3 

6 

4 

0 

0.7 

2FGLJ0309. 1+1027 


PKS 0306+102 

J03 12+01 

178.7 

-45.7 

0.8 

1 

0 

0 

1.8 

2FGL J0312. 6+0132 


PKS 0310+013 

J0320+41 

150.6 

-13.1 

0.3 

3 

3 

0 

0.8 

2FGL J0319. 8+4130 


NGC 1275 

J0329+22 

164.6 

-27.4 

0.8 

1 

0 

0 

1.8 

2FGLJ0326. 1+2226 


TXS 0322+222 

J0333-40 

245.2 

-54.1 

0.6 

1 

1 

0 

1.4 

2FGL J0334.2— 4008 


PKS 0332-403 

J0333+32 

158.4 

-18.9 

0.8 

1 

0 

0 

1.8 




J0341-02 

189.0 

-42.5 

0.8 

1 

0 

0 

1.8 

2FGL J0339.4— 0144 


PKS 0336-01 

J0342-25 

220.0 

-51.7 

0.6 

1 

1 

0 

1.4 




J0350+79 

130.5 

19.8 

0.6 

1 

1 

0 

1.4 

2FGLJ0354. 1+8010 


S5 0346+80 

J0351-21 

214.5 

-48.8 

0.8 

1 

0 

0 

1.8 

2FGL J0350.0— 2104 


PKS 0347-211 

J0402-36 

237.9 

-48.8 

0.2 

14 

7 

1 

0.6 

2FGL J0403.9— 3604 


PKS 0402-362 

J0422-01 

195.3 

-33.4 

0.6 

1 

1 

0 

1.4 

2FGL J0423.2— 0120 


PKS 0420-01 

J0427-60 

271.0 

-40.9 

0.6 

1 

1 

0 

1.4 

2FGL J0433.4— 6029 


PKS 0432-606 

J0428-38 

240.9 

-43.7 

0.2 

30 

11 

7 

0.5 

2FGL J0428.6— 3756 


PKS 0426-380 

J0441-00 

197.5 

-28.9 

0.4 

4 

2 

0 

1.0 

2FGL J0442.7— 0017 


PKS 0440-00 

J0448+1 1 

187.5 

-21.1 

0.4 

2 

2 

0 

1.0 

2FGL J0448.9+1 121 


PKS 0446+11 

J0456-23 

223.8 

-35.1 

0.2 

34 

11 

14 

0.5 

2FGL J0457.0— 2325 


PKS 0454-234 

J0503-01 

201.6 

-24.7 

0.6 

2 

1 

0 

1.4 

2FGL J0501.2— 0155 

4396 

PKS 0458-02 

J0504+04 

195.6 

-21.2 

0.3 

5 

4 

0 

0.7 

2FGL J0505. 5+0501 


PKS 0502+049 

J05 17+45 

163.1 

4.4 

0.6 

2 

1 

0 

1.4 

2FGL J0517. 0+4532 


4C +45.08 

J0520-36 

240.4 

-33.2 

0.6 

3 

1 

0 

1.4 

2FGL J0523.0— 3628 


PKS 0521-36 

J0525+16 

188.5 

-10.9 

0.8 

1 

0 

0 

1.8 




J0532-48 

255.0 

-32.6 

0.2 

9 

7 

0 

0.6 

2FGL J0532.0— 4826 


PMN J053 1-4827 

J0532+21 

184.7 

-6.3 

0.6 

3 

1 

0 

1.4 

2FGL J0534. 5+2201 


PSR J0534+2200 

J0533+07 

197.2 

-13.8 

0.3 

4 

3 

0 

0.8 

2FGL J0532. 7+0733 


OG 050 

J0536— 44 

250.2 

-31.6 

0.2 

16 

10 

1 

0.5 

2FGL J0538.8— 4405 


PKS 0537-441 

J0537-54 

262.8 

-32.4 

0.6 

1 

1 

0 

1.4 

2FGL J0540.4— 5415 


PKS 0539-543 

J0537+13 

191.9 

-9.3 

0.4 

5 

2 

0 

1.0 

Fermi J0539+1432 

3999 

TXS 0536+145 

J0539-34 

238.9 

-28.9 

0.8 

1 

0 

0 

1.8 

2FGL J0536.2— 3348 


1RXS J053629. 4— 334302 

J0604-70 

280.6 

-29.4 

0.6 

1 

1 

0 

1.4 

2FGLJ0601. 1-7037 


PKS 0601-70 

J0622+33 

179.7 

9.2 

0.2 

16 

13 

0 

0.5 

2FGL J0622. 9+3326 


B2 0619+33 

J0628-20 

228.8 

-13.9 

0.3 

3 

3 

0 

0.8 

2FGL J0629.3— 2001 


PKS 0627-199 

J0640+05 

206.5 

0.2 

0.6 

1 

1 

0 

1.4 

Fermi J0639+0548 

4224 

Nova Mon 2012 

J0646-30 

239.9 

-14.1 

0.8 

1 

0 

0 

1.8 

Fermi J0648— 3044 

3878 

PKS 0646-306 

J0652+45 

171.2 

18.9 

0.3 

4 

3 

0 

0.8 

2FGL J0654.2+4514 


B3 0650+453 

J0703-46 

257.1 

-17.5 

0.6 

1 

1 

0 

1.4 

2FGL J0701.7— 4630 


PKS 0700-465 

J0712+19 

197.6 

13.2 

0.4 

2 

2 

0 

1.0 

2FGL J0714. 0+1933 


MG2J07 1354+ 1934 

J07 17+71 

144.2 

27.7 

0.2 

14 

8 

3 

0.6 

2FGLJ072 1.9+7 120 


S5 0716+71 
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Table 1 

(Continued) 


Ackermann et al. 


1FAV ID 

/ 

(°) 

b 

n 

^68 

(°) 

N { 

< 

iy neg 

R* 

(°) 

LAT Assoc. 

ATel 

Assoc. 

J0720+33 

185.2 

20.1 

0.3 

4 

3 

0 

0.8 

2FGL J0719. 3+3306 


B2 0716+33 

J0724+14 

203.8 

13.6 

0.2 

10 

8 

i 

0.6 

2FGL J0725. 3+1426 


4C +14.23 

J0728-11 

227.6 

2.8 

0.4 

4 

0 

0 

0.9 

2FGL J0730.2— 1141 


PKS 0727-11 

J0730+36 

181.9 

23.3 

0.6 

1 

1 

0 

1.4 




J0740+54 

162.8 

28.8 

0.3 

11 

5 

0 

0.7 

2FGL J0742. 6+5442 


GB6 J0742+5444 

J0746+01 

217.9 

12.9 

0.8 

1 

0 

0 

1.8 

2FGL J0739. 2+0138 


PKS 0736+01 

J0747+24 

196.4 

22.6 

0.6 

1 

1 

0 

1.4 




J0759-56 

270.0 

-13.7 

0.6 

1 

1 

0 

1.4 

GRB080916C 



J0804+61 

155.3 

32.2 

0.6 

1 

1 

0 

1.4 

2FGL J0805. 5+6145 


TXS 0800+618 

J0806+52 

165.8 

32.4 

0.6 

2 

1 

0 

1.4 




J0808-08 

229.4 

13.0 

0.2 

6 

6 

0 

0.6 

2FGL J0808.2— 0750 


PKS 0805-07 

J08 11+02 

220.0 

18.9 

0.6 

1 

1 

0 

1.4 




J0831+04 

220.9 

24.1 

0.6 

1 

1 

0 

1.4 

2FGLJ083 1.9+0429 


PKS 0829+046 

J0839+00 

225.3 

24.4 

0.6 

1 

1 

0 

1.4 

2FGL J0839. 6+0059 


PKS 0837+012 

J0845+70 

143.6 

34.8 

0.6 

9 

1 

0 

1.4 

2FGLJ0841. 6+7052 


4C +71.07 

J0849-12 

238.6 

19.4 

0.3 

3 

3 

0 

0.8 

2FGL J0850.2— 1212 


BZQ J0850— 1213 

J0849+50 

168.3 

39.1 

0.2 

7 

6 

0 

0.6 

Fermi J0849+5108 

3452 

SBS 0846+513 

J0855+20 

206.4 

36.0 

0.3 

7 

4 

0 

0.7 

2FGL J0854. 8+2005 


OJ 287 

J0905-35 

259.9 

7.5 

0.6 

1 

1 

0 

1.4 

2FGL J0904.8— 3513 


NVSS J090442— 351423 

J0907-02 

232.7 

28.5 

0.8 

1 

0 

0 

1.8 

2FGL J0909.7— 0229 


PKS 0907-023 

J0908-50 

271.0 

-1.6 

0.6 

1 

1 

0 

1.4 

2FGL J09 10.4- 5050 


AT20G J091058— 504807 

J0908+60 

155.0 

40.0 

0.6 

1 

1 

0 

1.4 




J09 10+01 

228.9 

31.4 

0.6 

1 

1 

0 

1.4 

2FGL J0909. 1+0121 


PKS 0906+01 

J0920+44 

175.9 

44.8 

0.3 

12 

5 

0 

0.7 

2FGL J0920. 9+4441 


S4 0917+44 

J0924+28 

199.0 

44.5 

0.6 

1 

1 

0 

1.4 

2FGL J0924. 0+2819 


B2 0920+28 

J0948+01 

235.8 

39.1 

0.4 

3 

0 

0 

1.1 

2FGLJ0948. 8+0020 


PMN J0948+0022 

J0955+65 

145.7 

42.6 

0.6 

1 

1 

0 

1.4 

2FGL J0958. 6+6533 


S4 0954+65 

J0957-13 

251.7 

31.2 

0.6 

1 

1 

0 

1.4 

2FGL J0957.6— 1350 


PMN J0957- 1350 

J0957+24 

206.9 

50.9 

0.6 

1 

1 

0 

1.4 

2FGL J0956. 9+2516 


OK 290 

J1013+24 

208.4 

54.5 

0.6 

1 

1 

0 

1.4 

2FGL J1012. 6+2440 


MG2 J101241+2439 

J1015+04 

236.9 

46.8 

0.6 

1 

1 

0 

1.4 

2FGL J1016. 0+0513 


TXS 1013+054 

J1023-31 

269.6 

21.2 

0.6 

1 

1 

0 

1.4 




J1033+60 

147.8 

49.0 

0.2 

8 

6 

0 

0.6 

2FGL J1033. 9+6050 


S4 1030+61 

J1038-53 

283.6 

4.7 

0.3 

4 

3 

0 

0.8 

Fermi J1038-53 14 

3978 

PMN J1038— 5311 

J 1040+06 

240.9 

52.6 

0.6 

1 

1 

0 

1.4 

2FGL J1040. 7+0614 


4C +06.41 

J 1045+81 

128.4 

34.4 

0.6 

1 

1 

0 

1.4 

2FGL J1042. 6+8053 


S5 1039+81 

J 1046- 29 

272.6 

26.1 

0.6 

3 

1 

0 

1.4 

2FGL J1045.5— 2931 


PKS B1043— 291 

J1051+04 

245.7 

53.6 

0.6 

1 

1 

0 

1.4 




J1 100+00 

252.9 

52.7 

0.8 

1 

0 

0 

1.8 

2FGL J1058.4+0133 


4C +01.28 

J1 102+37 

180.8 

64.9 

0.6 

1 

1 

0 

1.4 

2FGLJ1 104.4+3812 


Mkn 42 1 

J1 110+34 

187.1 

67.1 

0.8 

1 

0 

0 

1.8 

2FGLJ1 112.4+3450 


TXS 1109+350 

J1 120-05 

265.7 

50.7 

0.3 

5 

3 

0 

0.8 

2FGLJ1 121.5-0554 


PKS 1118-05 

J1 124-64 

293.8 

-3.1 

0.8 

1 

0 

0 

1.8 

1FGLJ1 122.9-6415 


PMN J1 123-6417 

J1 129—19 

277.5 

39.5 

0.3 

6 

5 

0 

0.7 

2FGLJ1 126.6-1856 

3207 

PKS 1124-186 

J1 136—06 

271.8 

52.1 

0.8 

1 

0 

0 

1.8 




J1 145+39 

165.6 

71.3 

0.3 

5 

3 

0 

0.8 

2FGLJ1 146.9+4000 


S4 1144+40 

J1 152+49 

145.4 

64.8 

0.3 

8 

5 

0 

0.7 

2FGLJ1 153.2+4935 


OM 484 

J 1200+29 

197.9 

78.6 

0.2 

9 

6 

0 

0.6 

2FGLJ1 159.5+2914 


Ton 599 

J 1205+54 

136.6 

61.0 

0.6 

2 

1 

0 

1.4 

2FGLJ1208. 8+5441 


TXS 1206+549 

J1217+29 

190.3 

82.0 

0.6 

1 

1 

0 

1.4 

2FGLJ1217. 8+3006 


1ES 1215+303 

J 1224+21 

255.1 

81.6 

0.1 

111 

27 

82 

0.5 

2FGL J1224. 9+2122 


4C +21.35 

J 1229+02 

290.3 

64.5 

0.2 

34 

10 

14 

0.6 

2FGL J1229. 1+0202 


3C 273 

J 1240+04 

295.9 

67.1 

0.4 

4 

2 

0 

1.0 

2FGL J1239. 5+0443 


MG1 J12393 1+0443 

J 1246- 25 

301.5 

37.0 

0.2 

8 

6 

0 

0.6 

2FGL J 1246.7— 2546 


PKS 1244-255 

J1257-05 

305.5 

57.0 

0.2 

32 

10 

10 

0.5 

2FGLJ1256. 1-0547 


3C 279 

J 1300- 22 

305.6 

40.1 

0.6 

1 

1 

0 

1.4 

2FGL J1258.8— 2223 


PKS 1256-220 

J1303-64 

304.3 

-2.1 

0.3 

5 

0 

0 

0.9 

Fermi J1302-6350 

3085 

PSRB 1259-63 

J1312+55 

116.7 

61.6 

0.6 

1 

1 

0 

1.4 




J1313+32 

83.2 

82.9 

0.3 

3 

3 

0 

0.8 

2FGL J1310.6+3222 


OP 313 

J13 13+48 

113.4 

67.9 

0.2 

10 

7 

0 

0.6 

2FGL J1312. 8+4828 


GB 1310+487 

J1318-33 

309.5 

28.9 

0.6 

2 

1 

0 

1.4 

2FGL J1315.9— 3339 


PKS 1313-333 

J1329-70 

306.1 

-7.5 

0.4 

3 

2 

0 

1.0 

2FGL J1330.1— 7002 


PKS 1326-697 

J1331-55 

308.6 

6.6 

0.8 

1 

0 

0 

1.8 

2FGL J1329.2— 5608 


PMN J1329— 5608 

J1332-04 

321.7 

56.4 

0.2 

12 

7 

1 

0.6 

2FGL J1332.0— 0508 


PKS 1329-049 

J1333+04 

328.9 

65.4 

0.6 

1 

1 

0 

1.4 
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Table 1 

(Continued) 


Ackermann et al. 


1FAV ID 

; 

(°) 

b 

(°) 

^68 

(°) 

N f 


N f 

iV neg 

(°) 

LAT Assoc. 

ATel 

Assoc. 

J 1346+44 

94.7 

69.1 

0.2 

13 

8 

0 

0.6 

2FGL J1345.4+4453 


B3 1343+451 

J1348-29 

317.6 

31.6 

0.6 

1 

1 

0 

1.4 

2FGL J1351.3— 2909 


PKS 1348-289 

J1351-11 

325.6 

49.0 

0.3 

4 

4 

0 

0.7 

Fermi J1349- 1132 

3788 

PKS 1346-112 

J1353+30 

49.4 

76.1 

0.6 

1 

1 

0 

1.4 

2FGL J1350.8+3035 


B2 1348+30B 

J1419+35 

63.8 

69.3 

0.3 

4 

4 

0 

0.7 




J1428-41 

321.6 

17.5 

0.2 

14 

6 

4 

0.6 

2FGL J1428.0— 4206 


PKS B1424— 418 

J1459-35 

330.2 

20.4 

0.2 

6 

6 

0 

0.6 

2FGL J1457.4— 3540 


PKS 1454-354 

J1505+10 

11.5 

54.4 

0.1 

73 

26 

40 

0.3 

2FGL J1504.3+1029 


PKS 1502+106 

J1513-09 

351.4 

40.0 

0.1 

110 

19 

76 

0.5 

2FGL J1512.8— 0906 


PKS 1510-08 

J1517-31 

335.9 

21.8 

0.8 

1 

0 

0 

1.8 

2FGL J1513.6— 3233 


PKS 1510-324 

J1522+31 

50.0 

56.8 

0.2 

13 

7 

1 

0.6 

2FGL J1522. 1+3144 


B2 1520+31 

J1533-13 

352.2 

33.5 

0.2 

9 

9 

0 

0.5 

Fermi J1532— 1319 

3579 

TXS 1530-131 

J1554+13 

24.3 

45.2 

0.4 

3 

2 

0 

1.0 

2FGL J1553.5+1255 


PKS 1551 + 130 

J1626-25 

352.0 

16.0 

0.3 

6 

4 

0 

0.7 

2FGL J1625.7— 2526 


PKS 1622-253 

J1636+38 

61.3 

42.1 

0.2 

24 

6 

6 

0.6 

2FGL J 1635.2+38 10 


4C +38.41 

J1641+41 

65.0 

41.2 

0.6 

2 

1 

1 

1.4 

2FGL J1642.9+3949 


3C 345 

J1641+47 

73.8 

41.3 

0.6 

1 

1 

0 

1.4 

2FGL J1637.7+4714 


4C +47.44 

J1656+48 

74.6 

38.6 

0.6 

1 

1 

0 

1.4 

2FGL J1657.9+4809 


4C +48.41 

J1700+68 

99.7 

35.2 

0.3 

5 

3 

0 

0.8 

2FGLJ 1700.2+6831 


TXS 1700+685 

J1704-62 

328.1 

-12.7 

0.8 

1 

0 

0 

1.8 

2FGL J1703.2— 6217 


CGRaBS J1703— 6212 

J 1707+77 

109.3 

32.1 

0.8 

1 

0 

0 

1.8 




J1710+43 

68.5 

36.0 

0.2 

11 

10 

0 

0.5 

2FGL J1709.7+4319 


B3 1708+433 

J1715 — 33 

352.5 

2.8 

0.4 

2 

2 

0 

1.0 

2FGL J1717.7— 3342 


TXS 1714-336 

J1718-52 

337.5 

-8.5 

0.4 

2 

2 

0 

1.0 

Fermi J1717— 5156 

4023 

PMN J1717— 5155 

J1719+18 

39.8 

28.1 

0.4 

2 

2 

0 

1.0 

2FGL J1719.3+1744 


PKS 1717+177 

J1732-12 

12.1 

10.9 

0.3 

3 

3 

0 

0.8 

2FGLJ1733. 1-1307 


PKS 1730-13 

J1735+39 

64.3 

30.7 

0.3 

4 

4 

0 

0.7 

2FGL J1734.3+3858 


B2 1732+38A 

J 1740+27 

51.5 

26.9 

0.6 

1 

1 

0 

1.4 

GRB090902B 



J 1 74 1 +50 

77.1 

31.5 

0.6 

1 

1 

0 

1.4 

2FGL J1739.5+4955 


S4 1738+49 

J 1743+52 

79.5 

31.4 

0.4 

2 

2 

0 

1.0 

2FGLJ 1740.2+52 12 


4C +51.37 

J 1746+70 

100.7 

30.9 

0.3 

4 

3 

0 

0.8 

2FGL J1748. 8+7006 

3171 

S4 1749+70 

J1751+09 

34.6 

17.6 

0.6 

1 

1 

0 

1.4 

2FGL J1751.5+0938 


OT 081 

J1752-32 

357.3 

-3.3 

0.6 

1 

1 

0 

1.4 

Fermi J 1750- 3243 

4284 

Nova Sco 2012 

J1753+32 

58.0 

25.7 

0.4 

2 

2 

0 

1.0 

2FGL J1754.3+3212 


RX J1754. 1+3212 

J1759-48 

344.0 

-12.2 

0.6 

1 

1 

0 

1.4 

2FGL J1759.2— 4819 


PMN J1758— 4820 

J1803-39 

352.9 

-8.4 

0.2 

8 

6 

0 

0.6 

2FGL J1802.6— 3940 


PMN J 1802- 3940 

J1 807+78 

110.1 

28.7 

0.3 

7 

5 

0 

0.7 

2FGLJ 1800.5+7829 


S5 1803+784 

J1 825+56 

85.8 

25.9 

0.4 

2 

2 

0 

1.0 

2FGLJ 1824.0+5650 


4C +56.27 

J1827-52 

342.6 

-17.8 

0.6 

1 

1 

0 

1.4 

2FGLJ1825. 1-5231 


PKS 1821-525 

J1834-21 

12.3 

-5.9 

0.2 

10 

10 

0 

0.5 

2FGL J1833.6— 2104 


PKS 1830-211 

J1 850+32 

62.2 

14.3 

0.2 

7 

6 

0 

0.6 

2FGL J1848.5+3216 


B2 1846+32A 

J1 853+48 

78.3 

19.7 

0.6 

1 

1 

0 

1.4 

2FGLJ 1852.5+4856 


S4 1851+48 

J1 853+67 

97.6 

24.7 

0.2 

13 

11 

0 

0.5 

2FGL J1 849.4+6706 


S4 1849+67 

J 1909- 80 

313.9 

-27.5 

0.6 

1 

1 

0 

1.4 




J1910-20 

16.8 

-13.0 

0.5 

2 

0 

0 

1.3 

2FGLJ1911. 1-2005 


PKS B1908— 201 

J1914 — 35 

1.8 

-19.9 

0.6 

1 

1 

0 

1.4 

Fermi J19 13-3630 

2966 

PMN J 1913 — 3630 

J1924-21 

17.1 

-16.4 

0.3 

5 

3 

0 

0.8 

2FGL J1923.5— 2105 


TXS 1920-211 

J1956-38 

1.9 

-28.5 

0.6 

1 

1 

0 

1.4 

2FGL J1958.2— 3848 


PKS 1954-388 

J1956-42 

356.9 

-29.6 

0.6 

2 

1 

0 

1.4 

2FGL J1959. 1—4245 


PMN J1959— 4246 

J2001+44 

79.3 

7.3 

0.6 

1 

1 

0 

1.4 

2FGLJ2001. 1+4352 


MAGIC J200 1+435 

J2012+37 

74.6 

1.7 

0.6 

3 

1 

0 

1.4 

2FGL J2015.6+3709 


MG2 J201534+3710 

J2023+33 

73.0 

-2.1 

0.6 

1 

1 

0 

1.4 

2FGLJ2025. 1+3341 


B2 2023+33 

J2026-07 

37.0 

-24.6 

0.2 

15 

8 

0 

0.6 

2FGL J2025.6— 0736 


PKS 2023-07 

J2034+40 

80.0 

0.4 

0.4 

3 

0 

0 

1.1 

2FGLJ2032. 1+4049 


Cyg X— 3 

J2036+11 

55.7 

-17.3 

0.6 

1 

1 

0 

1.4 

2FGL J2035.4+1058 


PKS 2032+107 

J2056-47 

352.4 

-40.4 

0.3 

6 

4 

0 

0.7 

2FGL J2056.2— 4715 


PKS 2052-47 

J2 103+45 

87.1 

-0.7 

0.4 

2 

2 

0 

1.0 

2FGL J2102. 2+4546 


V407 Cyg 

J2 122-46 

353.8 

-44.9 

0.8 

1 

0 

0 

1.8 

2FGL J2125.0— 4632 

3808 

PKS 2123-463 

J2134-01 

52.3 

-36.5 

0.6 

1 

1 

0 

1.4 

2FGL J2133.8— 0154 

4333 

PKS 2131-021 

J2144+17 

72.2 

-26.2 

0.8 

1 

0 

0 

1.8 

2FGL J2143.5+1743 


OX 169 

J2148-75 

315.9 

-36.7 

0.3 

10 

5 

0 

0.7 

2FGL J2147.4— 7534 


PKS 2142-75 

J2154-30 

17.1 

-51.3 

0.3 

5 

0 

0 

0.9 

2FGL J2151.5— 3021 


PKS 2149-306 

J2155-83 

307.9 

-31.4 

0.6 

2 

1 

0 

1.4 

2FGL J2201. 9-8335 


PKS 2155-83 

J2 159+31 

85.1 

-18.5 

0.6 

4 

1 

0 

1.4 

2FGL J2157.4+3129 


B2 2155+31 

J2202+42 

92.5 

-10.4 

0.2 

15 

8 

1 

0.6 

2FGL J2202. 8+4216 


BL Lacertae 
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1FAV ID 

1 

(°) 

b 

n 

^68 

(°) 

N f 

< 

N f 

iV neg 

R s 

n 

LAT Assoc. 

ATel 

Assoc. 

J2202+50 

97.7 

-3.6 

0.3 

6 

4 

0 

0.7 

Fermi J2202+5045 

4182 

NRAO 676 

J2209-53 

339.8 

-50.2 

0.6 

1 

1 

0 

1.4 

2FGLJ2208. 1-5345 


PKS 2204-54 

J2214-26 

24.5 

-55.1 

0.6 

1 

1 

0 

1.4 

GRB090510 



J2229-07 

56.1 

-51.3 

0.5 

2 

0 

0 

1.3 

2FGL J2229.7— 0832 


PKS 2227-08 

J2233+1 1 

77.7 

-38.7 

0.3 

4 

3 

0 

0.8 

2FGL J2232.4+1 143 


CTA 102 

J2237-14 

48.3 

-56.3 

0.3 

5 

5 

0 

0.7 

2FGL J2236.5— 1431 


PKS 2233-148 

J2238-39 

0.4 

-60.0 

0.6 

1 

1 

0 

1.4 




J2247-06 

62.9 

-53.7 

0.8 

1 

0 

0 

1.8 




J2252-27 

24.5 

-63.6 

0.3 

7 

3 

0 

0.8 

2FGL J2250.8— 2808 


PMN J2250— 2806 

J2254+16 

86.2 

-38.2 

0.1 

168 

137 

123 

0.5 

2FGL J2253.9+1609 


3C 454.3 

J23 11+34 

100.5 

-24.1 

0.3 

9 

5 

0 

0.7 

2FGLJ231 1.0+3425 


B2 2308+34 

J2321+32 

101.7 

-27.1 

0.4 

2 

2 

0 

1.0 

2FGL J2322.2+3206 


B2 2319+31 

J2322-03 

76.8 

-58.6 

0.6 

1 

1 

0 

1.4 

2FGL J2323.6— 0316 


PKS 2320-035 

J2326+40 

105.8 

-19.6 

0.6 

1 

1 

0 

1.4 

2FGL J2325.3+3957 


B3 2322+396 

J2328-49 

332.4 

-62.3 

0.3 

27 

3 

2 

0.8 

2FGL J2329.2— 4956 


PKS 2326-502 

J2329-21 

45.3 

-70.5 

0.6 

1 

1 

0 

1.4 

2FGL J2330.9— 2144 


PMN J233 1—2148 

J2330+09 

91.8 

-48.5 

0.6 

3 

1 

0 

1.4 

2FGL J2327.5+0940 


PKS 2325+093 

J2332-66 

314.8 

-49.0 

0.6 

1 

1 

0 

1.4 

GRB090926A 



J2333-40 

348.4 

-69.1 

0.4 

2 

2 

0 

1.0 

2FGL J2336.3— 41 1 1 


PKS 2333-415 

J2345-15 

65.8 

-71.0 

0.1 

20 

17 

1 

0.4 

2FGL J2345.0— 1553 


PMN J2345- 1555 


Notes. The first column shows the FAVA identification number (ID), which is composed of the right ascension and declination of the source (J2000). 
The following columns show the Galactic coordinates and the statistical position error at 68% confidence R^. The systematic error on the source 
position is 0?1. Also shown are the total number of detected flares A' 1 , number of flares with detections at high energy A 7 :., , and number of flares 
corresponding to negative flux variations . The last columns show the associated Fermi-LAT source and the counterpart at longer wavelength found 
within a distance R s (see the text). For sources first announced via ATels we post the telegrams number. R s is derived as the 99% statistical error plus 
the systematic error. 
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Figure 6. Top panels show the relative variations of counts AA re i = (N — AI ex P)/A exp from the direction of Cyg X-3 (1FAV J2034+40) in weekly time bins, where 
N is the number of measured counts and IV exp is the number of expected counts from the average emission. Bottom panels show the significances of these variations. 
The left-hand plots are for photon energies above 100 MeV and the right-hand plots are for photon energies above 800 MeV. 

(The complete figure set (215 images) is available in the online journal.) 


the source position plus the systematic error r sys . R s was delib- 
erately chosen to be large, to include all possible counterparts. 
In cases for which more than one counterpart is found within 
R s , we consider the closest one. The values of R s for each source 
and the found counterparts are listed in Table 1 . We note that the 
associations were made purely on the basis of positional coinci- 
dence. We therefore caution that the associated sources should 
be considered as likely counterparts only. For a more confident 
source association, temporal and spectral information need to 
be considered. Additionally, the positional localizations of the 
FAVA sources could be improved by analyzing each source indi- 


vidually with standard likelihood techniques. This is beyond the 
scope of this paper. The associations were assigned as follows. 

1 . We searched for counterparts among the variable sources 
in the 2FGL catalog. We restricted the search to the 458 
2FGL sources that have a probability of less than 1% of 
being constant on monthly timescales. We find a variable 
2FGL source within the search radius for 170 of the FAVA 
sources. For those sources where no 2FGL counterpart was 
found, we searched for association with sources in the first 
Fermi source catalog ( 1FGL; Abdo et al. 2010a). The reason 
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is that sources that flared only once at the beginning of 
the mission might be detected in the 1FGL but not in the 
2FGL due to the increased integration time in the latter. We 
restricted the search to the 241 1FGL sources that have a 
probability of less than 1% of being constant, finding an 
association for one source. 

2. We searched for positional coincidences with Fermi- LAT- 
detected gamma-ray bursts 67 (GRBs). Even though GRBs 
have a typical duration from a few seconds to minutes, their 
emission is sometimes bright enough to be detected over a 
timescale of one week. We find a GRB within the search 
radius for four FAVA sources. These FAVA sources flared 
only once and we have verified that the flare occurred during 
the week of the GRB outburst. 

3. We searched for counterparts among LAT sources that were 
announced in Astronomer’s Telegrams 68 (ATels). These 
sources were found by the automated sky processing (ASP) 
used by the LAT Collaboration (Atwood et al. 2009). We 
found positional coincidences with 17 sources. 

We found LAT counterparts for 192 of the 215 FAVA sources. 
Most of the associated sources, 177, are AGNs. All associations 
found at higher Galactic latitudes (|7>| > 10) belong to this 
class. Among the AGN associations, 129 belong to the class 
of Flat Spectrum Radio Quasars (FSRQs) and 29 of them 
belong to the BL Lacertae (BL Lac) class. The number of 
gamma-ray-emitting BL Lacs is approximately the same as the 
number of FSRQs; FSRQs therefore flare intrinsically more 
frequently in the LAT energy range. This is in agreement with 
the observation that FSRQs are more variable in gamma-rays 
on monthly timescales (Ackermann et al. 2011b). Three of 
the FAVA sources are associated with non-blazar AGNs: two 
sources are associated with Narrow-line Seyfert galaxies, which 
were recently found to be variable gamma-ray sources (1FAV 
J0849+50, Donato & Perkins 2011; 1FAV J0948+01, Foschini 
et al. 2012), and one source is found coincident with the radio 
galaxy NGC 1275 (1FAV J0320+41; Kataoka et al. 2010). The 
remaining associated sources are AGNs of unknown type. 

4. GAMMA-RAY FLARES IN THE GALACTIC PLANE 

Of the 215 FAVA sources, 27 are detected at Galactic latitudes 
smaller than 10°; their positions are shown in the lower panels 
of Figure 5. We found associations to previously known LAT 
sources for all of them. The low-latitude FAVA sources can be 
grouped into a few categories. 

1. Sources associated with Galactic sources: four sources 
coincide with variable 2FGL sources which are associ- 
ated with Galactic sources. These are high-mass X-ray 
binaries Cyg X-3 (1FAV J308+41) andYsi +61 303 
(1FAV J0237+60), the Crab Nebula (1FAV J0532+21), 
and the nova V407 Cyg (1FAV J2103+45). In addition, 
three sources are found coincident with LAT sources an- 
nounced in ATels: the high-mass X-ray binaries PSRB 
1259-63 (1FAV J1303-64; Abdo et al. 2010f, 2011a) and 
two sources that are likely associated with novas (1FAV 
J1752-32, Cheung et al. 2012b; 1FAV J0640+05, Cheung 
et al. 2012c, 2012d). 

2. Sources associated with blazars: we find associations with 
variable 2FGL sources for 15 sources and to 1 1FGL 


67 http://fermi.gsfc.nasa.gov/ssc/observations/types/grbs/ 

68 https://www-glast.stanford.edu/cgi-bin/pub_rapid; 
http://www.asdc.asi.it/feratel (status 2012 November) 


source (1FAV J1 124—64), which were classified either as 
blazars or as AGNs of uncertain type in the respective 
catalogs. Additionally, three sources were found coincident 
with sources announced in ATels and associated with 
blazars (1FAV J0537+13, Orienti & D’Ammando 2012; 
1FAV J2202+50, Ciprini & Hays 2012; 1FAV J17 18-52, 
Chomiuk et al. 2013). 

3. Source with association of unknown type: one source is 
associated with counterparts of unknown type announced 
in ATels (1FAV J1038— 53; Ciprini et al. 2012). A compact 
radio source is found coincident with its position. 

Based on these associations, seven FAVA sources are located 
within the Milky Way. The source with the radio counterpart 
of unknown type might also be Galactic. Additionally, some 
of the associations made with blazars might turn out to be 
wrong. However, we can already infer statistically that most 
of the sources that we have not associated with Galactic sources 
are indeed extragalactic by calculating the number of expected 
extragalactic flares at low latitudes. The derivation relies on two 
assumptions. 

1. The majority of the sources at high Galactic latitudes are 
extragalactic. 

2. Extragalactic sources are isotropically distributed in the sky. 

We derive the number of extragalactic sources within 10° of 
the Galactic equator from the density of sources at latitudes 
greater than 30°. After considering the difference in solid angle, 
one expects 41 extragalactic sources at low latitudes. To take 
into account the reduced sensitivity for flare detection in the 
Galactic plane, we assigned random positions at low latitudes to 
the high-latitude sources. We determined the fraction of flares 
which would still be detected at the reduced sensitivity for 
each source. This results in an expectation of 24.6 variable 
extragalactic sources at Galactic latitudes smaller than 10°. The 
probability to detect 20 or fewer flares at low Galactic latitudes 
is 21%. The 20 sources not associated with Galactic sources are 
therefore compatible with being all extragalactic. No more than 
six of them can have a Galactic origin at >90% confidence. 

We note that two gamma-ray binaries LS 5039 and 1FGL 
J1018.6— 5856 (Abdo et al. 2009a; Hadasch et al. 2012; Ack- 
ermann et al. 2012c) are not detected by FAVA. Their orbital 
periods of 3.9 days and 16.6 days, respectively, result in average 
weekly flux variations below the sensitivity for flare detection 
by FAVA. The X-ray binary Cyg X-l is also not found by our 
analysis. The flare reported from this source by Sabatini et al. 
(2010) could not be confirmed by the LAT Collaboration. 69 

5. SUMMARY AND OUTLOOK 

We have presented the analysis tool FAVA, which searches 
for variable gamma-ray emitters in Fermi-LAI data. We used 
FAVA to search for gamma-ray flares on weekly timescales over 
the entire sky. From these flares, we derived a list of 215 flaring 
gamma-ray sources. A list of sources, their light curves, and their 
associations are available online. 70 We searched for positional 
coincidences of these sources with previously known LAT 
sources, finding counterparts for 192 sources. We associated 177 
sources with AGNs and find that FSRQs flare more frequently 
than BL Lacs. 


69 http://fermisky.blogspot.de/2010/03/lat-limit-on-cyg-x-l-during- 
reported.html 

70 https://www-glast.stanford.edu/pub_data/585 
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Twenty-seven of the FAVA sources are located at Galactic 
latitudes less than 10°. We associated seven of these to known 
Galactic sources. Among the remaining 20 sources, we found 
no evidence for new gamma-ray transients in our Galaxy. On 
the contrary, we showed that the majority of them are probably 
extragalactic. For 19 sources we find positional coincidence 
with AGNs. The one remaining source is associated with 
compact radio sources of unknown type. Future multi-waveband 
observations may reveal its nature. 

No flare was detected from a pulsar other than the Crab 
Nebula. It remains puzzling why this source is the only one 
of its kind to exhibit long-term variability and flaring behavior. 
We cannot confirm the hypothesis reported by Neronov et al. 
(2012) that gamma-ray variability might be common in young 
pulsars. We also detected no flares associated with previously 
undetected X-ray binary systems. These systems appear to 
be less efficient gamma-ray emitters than expected before the 
beginning of Fermi observations (Dubus 2007). We note that 
we detected no flares from the Galactic center region, which 
might have been expected if its gamma-ray emission was linked 
to accretion on the central black hole Sgr A* (Aharonian et al. 
2008). 

In the future we plan to apply FAVA on different timescales, 
and scan the gamma-ray sky for short-term flares on timescales 
of a few hours, and for long-term flux variations of a few months. 
Furthermore, we intend to run the analysis routinely to search 
for flares as soon as the LAT data are processed and sent to 
the Fermi Science Support Center. This will complement the 
ASP flare search currently used by the LAT Collaboration, and 
will help to alert the astrophysics community about gamma-ray 
flares in real time. 
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